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Abstract 

In this paper the effect of the magnetic field on the form of the 
equation of state and helium white dwarfs structure are discussed. 

The influence of the temperature and magnetic held on white dwarfs 
parameters have been investigated. The mass-radius relations for dif¬ 
ferent parameters were obtained. The occurrence of unstable branches 
in the mass-radius relation are presented for the temperature equals 
T = 3 10® K and for different values of the strength of magnetic held. 
Theoretical model of the star with two Landau levels is obtained. 

PACS numbers 64, 65, 97 

1 Introduction 

Properties of matter in strong magnetic field has been the subject of investi¬ 
gations in astrophysics of white dwarfs and neutron stars. It is motivated by 
the fact that magnetic fields of the order of 10® Tesla - 10® Tesla are known 
to exist in many cases of white dwarfs and neutron stars. White dwarfs are 
compact objects with masses comparable to that of the Sun, radii of the 
order of several thousands of kilometers and mean densities about 10® kg 
m“®. These stars no longer burn nuclear fuel. They are slowly cooling as 
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they radiate away their residual thermal energy. White dwarfs radii decrease 
with increasing mass. In accordance with theoretical predictions there is a 
critical valne of white dwarf mass known as Chadrasekhar mass which is 
predicted to be 


Mch = 5.8 Me: 


© 


where W is the fraction of electrons. This is the pressnre of degenerate elec¬ 
trons which snpported the star against collapse. The electrons are snpposed 
to be degenerate with arbitrary degree of relativity (in this paper the unit in 
which c = h = 1 was used) x = kF/me and it various with the density. He¬ 
lium ions provide the mass of the star and their contribntion to the pressnre 
is negligible ||l], |, i. 

The ions form a regnlar lattice which minimizes the energy. Detailed models 
of white dwarfs were given by Salpeter (1967) who derived conditions for 
their solidification and determined properties of the lattice. 

A nnmber of white dwarfs with strong magnetic fields was discovered (Kemp 
et al. 1970; Pntney 1995; Reimers et al. 1996) H, |] and extensively stnd- 
ied (Jordan 1992; Chanmngan 1992) [0, |]. 

In about 3-4 % of all white dwarfs magnetic fields have been detected rang¬ 
ing from 2 10“^ Tesla up to 10® Tesla. The magnetic white dwarfs where 
a surface magnetic field on the order of 10 Tesla ~ 10® Tesla and an inte¬ 
rior field of 10® Tesla ~ 10® Tesla are estimated. The magnetic field may 
canses considerable effects on the strnctnre of white dwarfs. Earlier work of 
Ostriker and Hartwick (1968) and recent calculations of Such and Mathews 
(2000) predicted an increase of white dwarf radii in the presence of inter¬ 


nal magnetic fields [|, 0]. In this paper the influence of the temperature 
and magnetic field on white dwarfs parameters have been investigated. The 
mass-radins relations for different parameters were obtained. 

This paper is organized as follows. 

In Sect.l are presented the general properties of white dwarfs. In Sect.2 the 
employed equation of state (EOS) is obtained in magnetic fields model with 
different temperatnres. Einally the EOS is nsed to determine the eqnilib- 
rinm confignrations of white dwarfs. The interesting fact is the existence of 
nnstable branches on the mass-radins relation for temperature T = 3 10® 
K and different values of the strength of magnetic field, which indicate that 
there are several possible evolutionary tracks. Einally, in Sect.3 the main 
implications of the resnlts are disenssed. 
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2 The white dwarfs in the magnetic fields model 

This paper presents the model of helium white dwarf in which the main con¬ 
tribution to the pressure comes from ultra-relativistic electrons. Corrections 
from finite temperature and magnetic field are also included. Having made 
the assumption that the ionised uniform helium plasma forms the interior of 
such an object one can say that the white dwarf matter consists of electrically 
neutral plasma which comprises charged ions {^He) and electrons. 

The Lagrangian density function in this model can be represented as the sum 

+ C-c + Cion + Cqed, (l) 

where £e, Cq describe the electron and gravitational terms, respectively. 
The Cion determines the ionized helium and Cqed is the Lagrangian density 
function of the QED theory. 

The electron part of the Lagrangian is given by 

Ce = “ rUelpei’e, 

where is the covariant derivative defining as — ieA^ and e is 

the electron charge. The Dirac equation for electrons obtained from the 
Lagrangian function C has the form 


{iYDi, - me)'ipe = 0 . 


The vector potential is defined as = {Aq = 0, Ai} where 


A* 


2 ^ilm^ B 


m 

0 • 


The gauge in which uniform magnetic field B lies along the z-axis was chosen 


B 


m 

0 


(0,0,5,). 


The magnetization is given by the contribution of the electrons and depends 
on the density of particles and antiparticles. Different from zero magne¬ 
tization generates internal molecular field 5, = fJ,M.z 7^ 0. The energy- 
momentum tensor can be calculated taking the quantum statistical average 

=< T^u >, 
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where 


T = 

fiiy 


2^ 

dgl^U 




In general in the presence of finite magnetization the pressnre is anisotropic. 
For the electrically charge particles, different eqnations of state for directions 
paraller and perpendicular to the magnetic field can be obtained. Thus the 
energy-momentum tensor has the form 


< T^u >= 


( e = p 

0 

0 

V 0 


0 

P 

0 

0 


0 0 \ 

0 0 

P 0 

0 P + pAizBz J 


( 2 ) 


This anisotropy in the pressure leads to a magnetostriction effect in the 
quantum magnetized gas of charged particles. In the classical case nonzero 
magnetization produces a flattening effect in white dwarfs and neutron stars 

Einstein equa- 


models ||TT| similarly like in the case of rotating stars 
tions (in isotropic case) leads to the standard Tolman-Oppenheimer-Volkov 
equations [T^. The equations describing masses and radii of white dwarfs 


are determined by the proper form of the equation of state. The aim of this 
paper is to calculate the equation of state for helium white dwarfs with the 
assumptions of finite temperature and in the presence of magnetic field. 

The properties of an electron in external magnetic field have been studied, 
for example by Landau and Lifshitz in 1938 [T^. In this paper the effects of 
high magnetic field on the equation of state of a relativistic, degenerate elec¬ 
tron gas is considered. The motion of free electrons in homogeneous magnetic 
field of the strength B perpendicular to the field is confined by the oscillatory 
force determined by the field B and is quantized into Landau levels with the 
energy nB, n = 0,1,.... In the case of nonrelativistic electrons the energy 
spectrum is given by the relation 


E 


n,Pz 


nuc + 


p 


2 

2 


2me 


where n = J + | + Sz, the cyclotron energy cOc = eBlrUe, Pz is the momentum 
along the magnetic field and can be treated as continuous. For extremly high 
magnetic field the cyclotron energy is comparably with the electron rest mass 
energy and this is the case when electrons became relativistic. Introducing 
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the idea of the critical magnetic field strength Be = rril/ \ e \ which equals 
4.414 X 10® Tesla, it is easy to distinguish between nonrelativistic and rel¬ 
ativistic cases, for B > the relativistic Dirac equation for electrons have 
to be used. The obtained dispersion relation now takes the form 

En,p, = \Jpl + ml + 2eB^n. (3) 

Along the field the motion is free, quasi one-dimensional with the modified 
density of states. The electron density of states in the absence of the magnetic 
field is replaced by the sum 



where the symbol 6no denotes the Kronecker delta [T^ and thus the spin 
degeneracy equals 1 for the ground {n = 0) Landau level and 2 for n > 1. 
The redefined density of states makes the distinctive difference between the 
magnetic and non-magnetic cases. The equation (||) implies that for n = 0, 
Eq = -|- ml whereas forn > 1 -f m*^. These relations indicate 

that the quantity m* defined as = rnl+2eBzn depends on the presence of 
the magnetic field and can be interpreted as an effective electron mass which 
is different from electron mass for n > 1. The number density of electrons at 
zero temperature is given by 

ne = S““^[2-5„o]e52pf. 

The maximum Landau level is calculated from the condition {p^Y — 0- 
can define the critical magnetic density ps, which denotes the limiting density 

Pb = 0.802 cm~^, (4) 

for densities lower than ps only the ground Landau level is present. 
Thermodynamic properties of a free electron gas in the magnetic fields at 
finite temperature have to be investigated. The temperature affects the elec¬ 
tron motion in external magnetic fields. Finite temperature and decreasing 
value of the strength of magnetic field tend to smear out Landau levels. 
The number density of electrons in the presence of a magnetic field can be 
expressed now as 

2iTi^'y 

Ue = - 5no\{hp,+ {z/t, 1 -k 2'yn) - Iofl-{z/t, 1 -k 27n)) (5) 
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where the Fermi integral 


i\,V± 


/■ (a + 

(m, tt j / / / \ 

J Ua+x^^u) , . 


( 6 ) 


^{Va+x‘^u^ ^ 

was used [|I^, [T^, 2 ; = fio/me, t = ksTo/me, u = zjt and 7 = BzjBc. 
Similarly to Fermi temperature which is given by 


Tp = Ep/kB = {me/kB)eF (for p < pb), 


one can define a magnetic temperature 


^ AEb me , , 

Tb = —r— = -j— ( \/1 + 2n; 


(^Y^l + 2nmax7 + 27 - + 2nmax7) 


(7) 


( 8 ) 


where is the energy difference between the np = Umax level and the 

= '^max + 1 level. For p = Pb these temperatures equal Tp = Tb. Knowing 
the value of magnetic temperature Tb one can describe how the properties of 
free electron gas changes at finite temperature and the presence of magnetic 
field. The influence of magnetic field is most significant for p < pB and 
T <Tb when electrons occupy the ground Landau level. In this case one can 
deal with the strong quantizing gas and magnetic field modifies all parameters 
of the electron gas. For example, for degenerate, nonrelativistic electrons the 
pressure is proportional to p^ and this form of the equation of state one can 
compare with that for the case B=0 for which P ~ When p > pB the 
Fermi temperature is still greater than magnetic temperature Tp electrons 
are degenerate and there are many Landau levels, now the level spacing 
exceeds fcpT. The properties of the electron gas are only slightly affected by 
magnetic field. With increasing temperature, there is the thermal broadening 
of Landau levels, when T > Tp, the free field results are recovered. For 
T ^ Tb there are many Landau levels and the thermal widths of the Landau 
levels are higher than the level spacing. The magnetic field does not affects 
the thermodynamic properties of the gas 


The total pressure of the system can be described as the sum of the 
pressure coming from electrons and ions plus small corrections coming from 
electromagnetic field 

P = Pe + Pion + PqED- 

The contributions from the same constituents form the energy density 


£ — £e p £ion + £QED- 
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The electron pressure and energy density are defined with the use of Fermi 
integral 


Pe = 


2'ym] 


~ ^no\iI-i, 2 ,+iz/t, 1 + 2'yn) + 1 + 2'yn)) (9) 


2^777,^ 

£e = ^^0 [2 ~ <^no](-fl,0,+ (^/^; 1 + 2^71) + 1 + 2'yn)). (lO) 

The ions are very heavy and treated as classical gas. The infiuence of mag¬ 
netic field for ions is very small so we have neglected this correction. The 
pressure and energy density dependence on the chemical potential po de¬ 
termines the form of the equation of state, which is calculated in the fiat 
Minkowski space-time. The obtained form of the equation of state is the 
base for calculating macroscopic properties of the star. In order to construct 
the mass-radius relation for given form of the equation of state the OTV 
equations have to be solved 


dP{r) 

dr 


-%(r) + P(r)) 


[m{r) -|- 47 rP(r)r^) 


( 1 - 


2Gm{r )' 
r ^ 


( 11 ) 


dm{r) 

dr 


47 rr^p(r). 


( 12 ) 


However, presence of strong gravitational field of the star causes the de¬ 
pendence of the temperature and chemical potential on the gravitational 
potential. In this paper together with the calculated form of the equation of 
state for completely ionised pure helium plasma the chosen values of central 
density pc changes from 10^ kg m~^ to 10^"^ kg m~^. For such parameters 
the value of magnetic field and temperature were limited to the following 
ranges T(0 — 10®) K, 7(0 — 1), respectively. Having solved OTV equations 
the pressure P(r), mass m(r) and density p(r) were constructed. To obtain 
the total radius R of the star the fulfillment of the condition P(r) = 0 is 
necessary. The results are presented in figures 1-10. Studying the properties 
of the star in the framework of the general-relativistic Thomas-Fermi model 
one can made the assumption that the temperature and chemical potential 
are metric dependent local quantities and the gravitational potential instead 
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of Poisson’s equation satisfies Einstain’s field equations. The metric is static, 
spherically symmetric and asymptotically fiat 


/ 0 0 0 \ 

0 -e^ 0 0 

0 0 —0 
V 0 0 0 —r‘^si'n?9 / 


(13) 


Its coefficients can be determined from Einstain equations and written as 
follows 

^ _ 2Gm(r) 
r 

The energy-momentum conservation 


TtM. = 0 


which for spherically symmetric metric ( [T^) can be written as 

du 2 dP 

dr P + p dr 


(14) 


together with the Gibbs-Duhem relation and with the assumption that the 
heat flow and diffusion vanishes [T^ give the condition 


u 

— = const. 
T 


(16) 


This implies that the temperature and chemical potential be come local met¬ 
ric functions 


T{r) = (16) 

p{r) = (17) 


where Tq and /tq are constants equal to the temperature and chemical po¬ 
tential at infinity. The temperature Tq may be chosen arbitrarily as the 
temperature of the heat-bath. Eirst equation in ([l6D is the well known Tol- 
man condition for thermal equilibrium in a gravitational field ||18|| . Eor given 
values po and Tq one can obtain self-consistency equations defining general- 
relativistic Thomas-Eermi equation. Now, the Eermi-Dirac distributions are 
defined as ||^ 


_ 1 _^_ 1 

(\/1+^^+27^^ —zo)/^0_j_^ (e^(^)/^\/l + x2+27n —2 q)/1q 


+ 1 


^(V l+a; 2 + 27 Ti —z(r))/t(T-)^^ 
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In the result instead of the OTV equilibrium equations together with the 
equation of state which was calculated in the flat space-time one can derived 
three self-consistency equations (|T^,|T^|T4D together with local form of the 


equation of state being now the function of r. This is of particular importance 
in the case of the strong gravitational flelds. 


3 Discussion 

In order to construct the mass-radius relation for white dwarfs the proper 
form of the equation of state have to be enumerated. For our needs we have 
chosen the pure helium plasma as the main constituent of the white dwarf 
interior. All calculations were performed at Unite temperature and different 
from zero magnetic held and compared with those of zero temperature and 
without magnetic held. In Figure 1 the applied form of the equation of 
state for magnetic and non-magnetic white dwarfs were presented. The zero 
temperature case is compared with the one obtained for temperature 6 x 10® 
K. For temperature equals zero there are clearly visible Landau levels which 
are smeared out when the temperature is different from zero and the strength 
of magnetic held is decreased. The density profile for helium magnetic and 
non-magnetic white dwarfs for different values of temperatures and central 
density pc equals 10^° kg was shown in Figure 2. As in the case of the 
EOS one can distinguish Landau levels for zero temperature case. The Fermi 
temperature for non-magnetic white dwarfs with the fixed value of the central 
density pc = 10^° kg is calculated and equals Tp = 1.023 x 10^° K. The 
estimated value of the magnetic held temperature for white dwarfs with the 
strength of magnetic held 7 = 1 (pc = 10^° kg m~^) is Tb = 1.27 x 10® 
K. Theoretical model of helium white dwarf with different electron masses 
which correspond to the ground Landau level (n = 0) and the one for n = 1 
is presented in Figure 2. The change of mass can be interpreted as the 
appearance two shells with distinctively different value of densities. For the 
presented form of the EOS the mass-radius relations were constructed. These 
relations are presented in Eigure 3. Of special interest is the case when the 
temperature equals 10® K because the appearance of unstable areas. The 
same is visible on the relation mass-central density (Eigure 4). There are 
marked points A, S, C on mass-radius and mass-central density relations. 
These points represents the stable configurations (A and C) whereas point 
B refers to the unstable one. Eigures 5 and 6 present sequences on the mass- 
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radius relations for fixed value of temperature. In Figure 5 the chosen value of 
the temperature equals T = 6 x 10® K whereas in Figure 6 the temperature 
equals T = 3 x 10® K . In both fignres the strength of magnetic field 
changes from 7 = 0 to 7 = 1. For comparison the zero temperatnre mass- 
radins relation obtained was presented. For increasing valne of the magnetic 
field white dwarfs increase their mass. However, the asymptotic value of 
Chadrasekhar mass is not exceeded. In the range of low densities together 
with increasing valne of the temperature the radii of these objects decrease. 
The mass-central density relations for helium white dwarfs are presented in 
Fignres 7 and 8 . Thick fnll cnrves stable branches, whereas thin fnll cnrves 
depict the unstable branches (Figure 8 ). The minimum and maximum of 
helium white dwarf star masses are represented by the extrema of M. These 
are point where the stable and dynamically nnstable branches merge. Stable 
branches are those with 4^ > 0. The existence of nnstable branches enables 

dpc 

different evolntion tracks. The temperatnre profile is presented in the Fignre 
9 for the fixed value of the central density pc- The temperature is changed 
slightly with the radius in agreement with the value of gravitational potential 
which is much smaller in comparison with neutron stars. White dwarfs with 
known values of masses and radii which collected in Table 1 are marked on 
the obtained mass-radins relations in Figure 10. In Fignre 10 the theoretical 
point m is marked. It represents the model of a star formed with one¬ 
dimensional electron gas. The changes of electron mass depend on the value 
of the Landau level n. The constructed configuration has distinctive shells 
with different valnes of the electron mass which is directly connected with 
the popnlation of Landan levels. For n = 0 (only the gronnd Landan level is 
popnlated) the electron mass is nnchanged whereas for u = 1 , one can deal 
with the modified by the magnetic field electron mass. 
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4 Figure captions. 

Figure 1 

The equation of state for different temperature and strength of magnetic field 
cases. 

Figure 2 

The density profile for helium magnetic white dwarf and for non-magnetic 
white dwarf for different temperatures and pc = 10 ^° kg/m~^. 

Figure 3 

The mass-radius relation for white dwarf with different temperature cases. 
Figure 4 

The mass M dependence on the central density pc for different temperature 
values. 

Figure 5 

The mass-radius diagram for helium white dwarfs for temperature T = 6 10® 
K and the strength of magnetic field 7 = 0.1 — 1. The solid line denotes the 
mass-radius relation for T = OiF. 

Figure 6 

The mass-radius diagram for helium white dwarfs for temperature T = 3 10® 
K and the strength of magnetic field 7 = 0.1 — 1. 

Figure 7 

The mass M dependence on the central density pc for temperatures T = OiF 
and T = 6 10® K. 

Figure 8 

The mass M dependence on the central density Pc for temperatures T = OK 
and T = 3 10® K. 

Figure 9 

The temperature profile for different central density. 
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Figure 10 

The mass-radius relation for white dwarf with different values of temperature 
with marked observational white dwarfs. The point m denotes the theoretical 
model of the star with Landan levels nL = 0,nL = l which are clearly visible 
in the Fignre 2. The mass and radins for this confignration is calcnlated 
R = 0.92 (lOOi?/i?0), M = 1.18 Mq. 
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Object 

M[Mq] 

p mu 

K p 

L /t(7) J 

Ref. 

SyriuszB 

1.053 ±0.028 

0.74 ±0.06 




40 Eri B 

0.48 ±0.02 

1.24 ±0.05 


1 


Stein 2051 

0.50 ±0.05 or 0.72 ± 0.08 

1.15 ±0.12 


T 


DA 

0.45 ±0.015 

1.85 ±0.08 


2C 


GD448 

0.41 ±0.01 

1.8 


21 


GD191-B2B 

0.51 

2.0 


22 


Feige 

0.50 

1.8 


22 


WD2218+706 

0.40 

3.3 


52 


REJ0558+165 

0.58 

1.9 


22 


REJ1738+665 

0.53 

2.4 


22 



Table 1: The masses and radii for some typical white dwarfs. 
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